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The mobility of damaged chromatin regions in the nucleus may affect the probabUity of mis-repair. In this 
work, hve-cell observation and distance tracking of GFP-tagged DNA damage response protein MDCl was 
used to study the random-walk behaviour of chromatin domains containing radiation-induced DNA 
double-strand breaks (DSB). Our measurements indicate a subdiffusion-type random walk process with 
similar time dependence for isolated and clustered DSBs that were induced by 20 MeV proton or 43 MeV 
carbon ion micro-irradiation. As compared to normal diffusion, subdiffusion enhances the probability that 
both ends of a DSB meet, thus promoting high efficiency DNA repair. It also limits their probability of 
long-range movements and thus lowers the probability of mis-rejoining and chromosome aberrations. 

DSB are a major threat for genomic integrity. However, some pathways for DSB repair can mis-rejoin DNA 
ends from different break sites, thus creating chromosome aberrations. The influence of initial distance 
between DSB sites on the probability of mis-rejoining is under debate, as well as the possibility of 
mechanisms bringing DSB sites into closer vicinity during damage processing'. Estimates on break interaction 
distances based on biophysical modelling vary between 100 nm and several \im^^. Simulation models for DNA 
repair based on normal diffusion of ends tend to overestimate the fraction of unrejoined DSB after long repair 
times^. 

Experimental data on the mobility of break sites in mammalian chromosomes are inconclusive. Mostly, they 
rely on labelling the chromatin region containing the DSB by visualizing specific chromatin modifications known 
to surround the breaks site, such as y-H2AX, or proteins known to bind chromatin in the vicinity of break sites, 
such as MDCl. The labelled chromatin regions microscopically appear as ionizing radiation-induced foci (IRIF). 
Analyzing fixed cells, inferences on mobility are made by investigating the development of foci number and size 
over time, or their position with regard to the sites of damage induction'' In live-cell analysis, also the behaviour 
of individual foci over time can be observed and several authors measured trajectories in relation to the origin 
(single particle tracking" '"'). In this case, thorough correction for translational and rotational movement of the 
whole nucleus, as well as its deformation, is paramount, since otherwise the measured mean squared displace- 
ments MSD(/lt) = <r^(At)> reflect movement of the nucleus rather than that of the foci within the nucleus. 

A plot of MSD versus observation time allows determining if the particle observed moves via normal diffusion, 
in which case the MSD is linearly related to time and a diffusion coefficient D can be determined'"'"'. Assuming a 
linear dependence of MSD on time, recently a large diffusion coefficient was determined" which was difficult to 
reconcile with the long-term stability of microscopic foci patterns after irradiation with ions in geometric 
patterns'. In addition, others reported for regions of undamaged chromatin that mobility is best described by 
subdiffusion (reviewed by''). In subdiffusion, MSD does not follow linear time dependence, but obeys 



MSD{At) = {r^{At))-- 



r{i+c 



-Af 



(1) 



with Drj_ being a generalized diffusion coefficient and oc being the anomality parameter quantifying the deviation 
from normal diffusion. In the case of subdiffusion, a < 1, while in normal diffusion, a = l.Y represents the 
Gamma function and d is the dimension of the analyzed space. 

To characterize the mobility of damaged chromatin regions, we chose to investigate distance changes between 
two IRIF over time rather than the MSD, using the standard deviation of the distance changes a^(Al). While 
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measuring distance changes has been introduced by others'" " due to 
the obvious advantage over single-particle tracking in case of any 
translational and rotational movements of the nucleus, our method 
using the (f(Al) is even largely invariant in the case of additional 
changes of nucleus size. 

Results 

IRIF mobility follows a subdiffusion process. U20S cells express- 
ing GFP-tagged repair protein MDCl were irradiated in a matrfac 
pattern with single carbon C^C) ions per point. The distances Z, of 
adjacent MDCl foci pairs numbered by i were monitored for several 
hours post-irradiation (Fig. 1). The distance changes Ali(At) during 
different time intervals /It (i.e. Ali(At) = li(tj + At) — li(tj) for all time 
points tj of the time series) were determined. 

For each cell, we determined the standard deviation a^(Al(At)) of 
the changes in foci distances Alj(At) between neighbouring MDCl 
foci as a measure for the underlying random walk process (see SI). 
Then the time-ensemble average of the "neighbouring" a^(At) from 
all analyzed carbon irradiated cells was calculated for each At and 
plotted versus AtinsL double-logarithmic plot (Fig. 2, filled squares). 
A linear increase of log(o"^) with \og{At) is evident, which demon- 
strates a power-law dependence of on At (c£ Supplementary 
equation (SI)): 

H{'^') = ^og(^-I^^^+c,logiAt) (2) 

We obtain a = 0.49 ± 0.05 (SEM) and D0.49 = (3.7 ± 0.5) X 
10"* iim^/s" *' for the next neighbour foci {I ~ 5 |J,m) of the carbon 
irradiated cells (R^ = 0.97, reduced %^ = 1.6). The value obtained for 
a significantly differs from a = 1 demonstrating that the mobility of 
the IRIFs follows a subdiffusion process^". Our data do not support 
models based on constrained diffusion^', since we did not see indica- 
tions for a limit to the radius within which the IRIF can diffuse, in 
spite of the long time range (up to IC s) observed. In addition, the 
best least-square fits of constrained diffusion" (R^ = 0.91, reduced 
= 5.0) or a combination of constrained diffusion and normal dif- 
fusion of the confinement region (R^ = 0.95, reduced = 2.8) did 
not fit the experimental data (Supplementary Fig. S2). 

To examine the influence of the initial distances on the determined 
dynamics, the cells were additionally analyzed using the distances 
between second next foci neighbours (/ = 10 |xm). As shown in Fig. 2 
(open circles), the analyzed from both data sets are very similar for 
At < 180 s, but increase faster in the case of second next neighbours 
at /ft > 180 s, hinting at non-corrected additional relative foci move- 
ments which we attribute to gross morphological deformations of 
cell nuclei. We determined a = 0.58 ± 0.05 (SEM) andDo.ss = (3.1 ± 
0.5) X 10-* [irnVs""'. 



■ C - next neighbours 
0 C - second next neighbours 




l,(tj) . l,(tj+At) 

* 5 |jmi * 5 |jm 



Figure 1 | U20S cell irradiated witli carbon ions in a 5 x 5 fini' matrix 
pattern, visible under the fluorescence microscope by GFP-tagging of 
DNA repair protein MDCl which accumulates at the radiation induced 
double-strand breaks ("foci"). The distance Z, between neighbouring foci 
is recorded at several time points tj for analyzing the distance changes A Z, 
over different time intervals A t. 
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Figure 2 | Double-logarithmic plot of the squared standard deviations tr^ 
(±SEM) of the distance changes Jl(Jt) between neighbouring MDCl 
foci (fUled squares) and second-next neighbours (open circles) in the 
nuclei of cells irradiated with carbon ions (657 foci pairs from 58 cells, 30- 
60 time frames per sample). The data are fitted with the power-law 
function (equation (2)). 



Comparison between time-series recorded early (3-10 min) and 
late (80-100 min) after irradiation shows that the subdiffusion of the 
MDCl foci is time-invariant, justifying the used method of 
ensemble- and time-averaging as it can be done for any ergodic 
process (Supplementary Fig. S3). 

The anomality of the diffusion of MDCl foci can also be recog- 
nized in the histogram of distance changes /)Z, for a single time 
interval At obtained for the carbon-irradiated samples (Fig. 3). 
Normal diffusion would lead to a Gaussian distribution of the prob- 
ability density ^(/fZ, At) of the Ali(At) values for initial distances 
much larger than the diffusional length^". The measured distribution 
function differs significantly from a Gaussian distribution, but is well 
fitted by the distribution function for subdiffusion (adapted from^") 
when using an anomalous diffusion exponent a = 0.5 (see Fig. 3, 
where zlt = 60 s) and applying a two dimensional random walk 
situation for distance analysis and thus substituting Al and 

2dD^. 

p,.Mo^^g: .._.'-."' (3) 



\/8^A^ fr^ /c!r ( 1 - ^ (Zc + 1) 



Vv/2dD^zlt« 



The differences to the Gaussian are visible both at very small dis- 
tances, where a cusp-shaped histogram is obtained as predicted by 
the subdiffusion model (see plot on linear scale of Fig. 3fl), as well as 
in the tails, where the distributions are significantly higher than 
predicted by a best Gaussian fit (log-plot of Fig. 3fc). We conclude 
from our live-cell imaging and distance tracking over time periods 
between a few seconds and several hours after irradiation that IRIF 
do not diffuse normally but follow a subdiffusion process. 

LET dependence of apparent foci mobility. To investigate whether 
IRIF diffusion depends on the density of damage sites along the 
irradiation track, the cells were also irradiated in a matrix of Ax = 
Ay = 6 \xm with 32 protons (20 MeV at cell position, linear energy 
transfer LET = 2.6 keV/|.im) per point. The number of protons was 
chosen to create on average 3 DSB per matrix point so that within the 
depth of view each detected IRIF contained on average only one DSB, 
while after carbon irradiation several DSB are expected to occur in 
each IRIF. As in the case of carbon-irradiated cells, a linear increase 
of log(t7^) with log(/ff) is evident (Fig. 4), demonstrating a power-law 
dependence of on A t. For comparison, also the data obtained after 
carbon ion irradiation are shown. 
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Figure 3 | Histogram of distance changes {Jl{Jt = 60 s) between two neighbouring foci after carbon irradiation, fitted with the distribution 
function for normal and anomalous diffusion on linear (A) and logarithmic scale (B) (bin size = 0.03 fim, number of data points Ntatai = 5656). 
SE-error bars are taken as \/N given through Poisson statistics. 



We determined oi = 0.50 ± 0.05 (SEM) andDo.50 = (8.1 ± 1.8) X 
10"* [.im^/s"-^" for the proton-irradiated cells, showing that the time 
dependence of foci after carbon and proton irradiation is very sim- 
ilar. As the anomality parameter a for proton and carbon irradiation 
(a = 0.50 ± 0.05, and a = 0.49 ± 0.05, respectively) is the same 
within the errors, the number of induced double-strand breaks 
within one focus does not seem to influence the degree of anomality 
of the foci diffusion. The absolute values, however, differ by a 
factor 2.2. This difference may be attributed to the fact that the 
dynamics observed in proton irradiated cells is mainly that of an 
isolated, single DSB (number of DSB per IRIF Np ~ 1), whereas 
due to clustered DSB induction and poor axial microscopic resolu- 
tion a number of Nc > 1 foci contribute to each detectable IRIF after 
carbon irradiation. Thus, the diffusional behaviour observed is that 
of the centre of mass of all DSB in this IRIF, where contributions of 
individual DSB can neutralize each other. In a pure statistical 
approach we expect a^pla^c — Nc/Np (see Supplementary Informa- 
tion "Random walk of the centre of mass of N DSB"). The measured 
ratio of 2.2 is approximately the number of DSB per micrometer as 
calculated in" multiplied with the approximate depth of field of the 
used epifluorescence microscope. 
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Figure 4 | Double-logarithmic plot of the squared standard deviations 
(±SEM) of the distance changes Jl(Jt) between neighbouring MDCl 
foci in the nuclei of cells irradiated with carbon ions (filled squares, 657 
foci pairs from 58 cells) and protons (open circles, 32 p per irradiation 
point, 99 foci pairs from 12 cells), each 30-60 time frames analyzed per 
sample. The data are fitted with the power-law function (eq. (2)). 



To experimentally verify this assumption, we rationalized that 
increasing the microscopic resolution would allow to differentiate 
between the IRIF of individual DSB. With smaU-angle irradiation it is 
possible to visualize with good lateral resolution how ion-induced 
IRIF form along the track''. Due to technical limitation, we were not 
able to perform live-cell analysis after small-angle carbon irradiation. 
We decided, therefore, to reinvestigate data obtained in earlier work 
(Du et al. 2011") by NBSl foci detection in carbon-irradiated cells 
fixed at various post-irradiation intervals. To allow a comparison of 
these data to our live-cell data, we introduce the "apparent diffusion 
coefficient" Dapp, which is defined as the ratio of the mean squared 
displacement MSD(At) (equation (1)) and the corresponding time 
At, corrected for the dimensionality d of the random walk: 

MSD{At) _ 2D„ 
2d At " r(l + a 



DappiAt) 



-Af 



(4) 



Our live-cell analysis has d = 2 due to the two-dimensional analysis 
of the centre of foci in the focal plane of the epifluorescence micro- 
scope, while the analysis in" has been performed in one dimension. 
The apparent diffusion coefficient D^pp is equal to the average dif- 
fusion constant that would result in the same MSD if normal dif- 
fusion had taken place. It is evident (Fig. 5) that the mobility of IRIF 
induced by small-angle carbon irradiation compares well to the 
mobility determined in the present work for proton-induced IRIF. 

Discussion 

In the present work, we demonstrate that radiation-induced MDCl 
foci, which are assumed to reflect the chromatin region surrounding 
one DSB or several clustered DSB, exhibit a subdiffusive mobility in 
the time range 5 s £ z)t £ 10* s, with a = 0.5. Subdiffusive behaviour 
of chromatin regions has been described in yeast^^, E. coli chromo- 
somes^''^* and, at least in the time scale of up to 100 s, telomeric 
regions in mammalian U20S ceUs^^'^*". Others, however, reported 
constrained diffusion of genomic loci (reviewed by^'). Recenfly, 
Albert et al. (2012) proposed that genomic loci in general undergo 
subdiffusion at all time scales explored, and that this effect may be 
overlooked when plotting MSD on a linear scale, where constrained 
and subdiffusive mobility lead to very similar curves". 

The underlying mechanisms leading to subdiffusion have mainly 
been investigated in the context of anomalous mobility of macro- 
molecules such as proteins and microinjected beads in nucleoplasm 
and cytoplasm (e.g.^^'^") where subdiffusion appears to be a con- 
sequence of crowding-induced viscoelasticity'""'^'. Modelling results 
suggest that subdiffusive motion of chromosomal loci in E. coli is a 
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Figure 5 | The apparent diffusion coefficient Dapp(At) -- 



MSD(At) 



ofDSB 



induced by different ionizing radiation (A) and in comparison with that 
of centromeric, telomeric and interstitial chromosomal loci (B) (with 
standard errors). C and p are the data sets from this work (cf. Fig. 4). 
(A) "C (Du)" is the lateral diffusion data of NBSl-foci in BJl-hTERT cells 
irradiated with C ions (~40 MeV) at SNAKE" and "Ni (Jakob)" and "y 
(Krawczyk)" are the diffusion data of GFP -tagged 53BPl-foci in U20S cells 
irradiated with 6 MeV Ni ions'' and with y", respectively. (B) Telomere 
data extracted from Bronstein et aP^., Jegou et al"., De Vos et al". and 
Krawczyk et aP'. Apparent diffusion coefficient of centromeric and 
interstitial chromosomal loci calculated from Chubb et aP*. and Krawczyk 
et aP'. Dapp is calculated using eq. (4) with corresponding dimension d and 
calculating MSD according to eq. (1) in case of distance measurements. 

consequence of both viscoelasticity of the nucleoplasm and the 
restrictions on monomer mobility caused by neighbouring mono- 
mers in a polymer^'. Also transient immobilization by binding to 
obstacles may contribute to subdiffusive behaviour™"'". 

The most pertinent characteristic of subdiffiision is that apparent 
mobility becomes slower at longer observation times. To detect sub- 
diffusive behaviour it is therefore important to avoid overestimation of 
mobility at long observation intervals which may be caused by insuf- 
ficient correction for mobility of the reference system, i.e. the nucleus 
in the case of chromosomal loci. We derived D^pp from a variety of 
data on mobiUty measurements of chromosomal loci available in the 
literature (Fig. 5a + b.). We note that those few mobiUty measure- 
ments which did not show a reduction of D^pp with observation time 
were performed with single particle tracking (Fig. 5a "Ni (Jakob)"", 
Fig. 5b "telomere (Bronstein)"^'' at long /it). In contrast, distance 
tracking as used in our work, intrinsically corrects for translational 
and rotational movements and, when using the standard deviation of 



distance changes as it is used here, even largely corrects for a change of 
the nucleus size. It is important to note that reduced apparent mobility 
at longer observation periods is not a consequence of altered condi- 
tions with increasing post- irradiation time (e.g. due to repair, or 
unphysiological conditions) since the mobility at short observation 
intervals is similar whether investigated within the first minutes after 
irradiation or after one hour (see Supplementary Fig. S3). 

What are the consequences of subdiffusive mobility? As discussed 
by Weber (2010), "subdiffusive molecules take longer than expected 
to reach distant targets but they are also more likely to retrace local 
regions of space than a freely diffusing particle"'''. Thus, the encoun- 
tered subdiffusion predicts long-range movements of chromosomal 
loci and DSB to be rare, even at long time intervals. Figure 5a shows 
that our measurements for MDCl and NBSl foci" mobility compare 
well with data on 53BP 1 foci determined recently after y irradiation"". 
Since contradictory data were published on whether IRIF can travel 
Jim distances and coalesce to form repair centres''"'"''^, it is inter- 
esting to compare IRIF mobility to that of undamaged chromatin 
regions. As evident from Fig. 5b, these data agree also well with a 
large set of published data on mobility of telomeric, centromeric, or 
interstitial chromosomal loci. In spite of the use of different methods 
and experimental systems, all available mobility data are surprisingly 
similar. Whether, as suggested by Krawczyk et aF'., IRIF and telo- 
meric loci exhibit generally a slightly enhanced mobility as compared 
to interstitial loci remains to be investigated in detail, also taking into 
account differences in foci substructures. 

In any case, we conclude that any model invoking movements of 
foci over several |tm within a few hours cannot be explained by 
typical mobility of chromosomal loci. Since the average distance of 
travel is small, the probability of joining wrong ends from DSB that 
initially were separated by several micrometers is low. For example, 
with typical values for D^pp smaller than 10"^ |tm^/s at zit = 3 h 
(Fig. 5b), the average sampled volume has a radius of y/ldDappAt 
= 770 nm {d = 3) which is less than 0.5% of the total volume of the 
cell nucleus of about 700 [iml 

Our study investigated the mobility of IRIF surrounding break 
sites rather than the mobility of individual DNA ends. Little is known 
on the extent of tethering of both ends of a DSB site, resulting either 
from the original chromatin structure or from fast DSB detection and 
response mechanisms'*^. There must be a certain probability of sepa- 
ration between ends, since otherwise mis-rejoining and generation of 
aberrations would not occur. For our considerations, we assume that 
the single ends underlie subdiffusion of similar speed as the DSB site 
as a whole. If there was independent movement of both single ends at 
considerably faster velocity, the monitored IRIF would be expected to 
frequently dissociate into two foci, which was not seen, although rare 
dissociation events cannot be excluded from the limited statistics of 
the experiment. Furthermore, the DNA ends are still connected to 
the chromatin, which is found to foUow subdiffusion also in the case 
of telomeric ends, with a mobility that agrees well with the IRIF 
mobility (cf. Fig. 5b). Computer simulations demonstrated that sub- 
diffusion enhances the probability of finding nearby targets as com- 
pared to normal diffusion for initial separations R smaller than the 
so-called crossover radius R^, thus facilitating cellular processes such 
as protein complex formation with high probability"'. A target size a 
= 2 nm to find a broken DNA end approximated by the diameter of 
the size of the DNA end, a nucleoplasmic viscosity of 2 X 10"'' Pa"s 
at 37°C and an anomalous diffusion coefficient a = 0.5 lead to a 
search time-dependent crossover radius R^ [jim] ~ 0.5 X tsearcii"* 
(with tjearch in [s]) of more than 5 jim for search times of 1-3 hours 
(calculated according to""). The radius may be even enhanced by a 
factor 2"" = 1.2 because both ends of a DSB are subject to subdiffu- 
sion. Provided that the distance of the two ends of a DSB that had 
separated from each other after damage induction is smaller than R^, 
the probability that the ends encounter again would thus be larger in 
the case of subdiffusion than for normal diffusion. This means that 
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subdiffusion supports the rejoining of ends during repair of DNA 
double-strand breaks. By labelling both sides of an enzyme-induced 
break site with different fluorescent molecules, Soutoglou et al. 2007 
demonstrated typical separation distances in the range of = 
100 nm^^, i.e. much smaller than the crossover radius R^. of more 
than 5 |im estimated here for several hours of search/repair time. 
Thus, the rejoining probability (within a few hours after irradiation) 
can be estimated to be larger than 90% for subdiffusion (using the 
model oP^), while normal diffusion would only rejoin broken ends 
with initial separation R ~ 100 nm with maximum 2% probability 
iP(R) = a/R'^). 

When there are more than one DSB within the initial separation 
distance R of the two ends of a single DSB, the probability to connect 
two ends would be the same for all DNA ends, thus resulting in a 
large fraction of translocations. Note that the ratio of correctly 
rejoined ends and translocations would be the same under subdiffu- 
sion or normal diffusion random walk, but many more single ends 
would be left in case of a normal diffusion. In the case that two DSB 
are created at a distance larger than R, however, the ratio of correctly 
rejoined ends to translocations increases more for subdiffusion than 
for normal diffusion. 

In conclusion, we propose that the subdiffusional behaviour of 
DSB sites guarantees a high probability of rejoining the ends of an 
individual DSB while reducing the probability of chromosomal aber- 
rations produced by rejoining the ends of different DSB or of leaving 
a large fraction of DNA ends unrepaired. 

Methods 

Cell culture and irradiation. Generation of a stable transfectant (clone Fl) of human 
osteosarcoma cells {U20S) with pEGFP-MDCP'* was seeded into custom-made cell 
containers'^ in HEPES-buffered, phenol red-free medium supplemented with 

0. 25 mM Trolox at 37 C and 5% CO2 (72 hours before irradiation). 

The irradiation was performed at the microbeam facility SNAKE^^'^^ at the Munich 
14 MV tandem accelerator with one carbon ion (43 Me V, LET ^ 370 keV/|im) or 32 
protons (20 MeV, LET — 2.6 keV/|im) per point in a 5 X 5 |j.m^ matrix, with an 
accuracy of 0.58 ]im in x and 0.86 |j.m in y'^ (see Fig. 1). One single 43 MeV carbon 
ion results in approximately 20 DSB along its approximately 7 |j,m track through the 
nucleus^, and the lowest number of protons necessary to obtain observable foci, which 
was about 30 protons per point, means an average number of three DSB that are 
distributed randomly along the ion track across the cell nucleus that has a height of 
about (7.2 ± 0.3) \xm. During irradiation and the subsequent observation the cells 
were covered by cell culture medium and the temperature was kept constant at 37''C. 

Distance tracking with 2D Fluorescence Microscopy and quantitative dynamics 
analysis. Live cell observation started approximately 200-400 s after irradiation and 
time series of different interval length At (e.g. 5 s, 1 or 10 minutes) and total length 
ttotai (20 minutes up to several hours) were recorded. Motion analysis was performed 
with the open source image analysis software Image} (http://rsb.info.nih.gov/ij) using 
the SpotTracker plugin (http://bigwww.ep.ch/sage/soft/spottracker/), developed by 
Daniel Sage et al'^. for tracking fluorescent particles in dynamic image sequences with 
subpixel resolution. After visually verifying that the foci were tracked correctly, the 
two-dimensional distance l(t) of all neighboring foci pairs (as well as second next 
neighbors) were calculated from the x-y coordinates gathered by the SpotTracker 
algorithm and the standard deviation a of the distance changes Alj was determined 
for different time intervals At for analysis of the foci dynamics. For more details see 
Supplementary Information. 
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